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I. INTRODUCTION 


The marine boundary layer is that region in the first 
kilometer over the sea which is defined principally by tur- 
bulent transfer of momentum and heat. In the past decade 
there have been concentrated efforts to make measurements 
in this region and to define physical expressions describing 
the turbulent intensities. The marine environment presents 
a number of problems, especially in accessibility, platform 
stability, and sensor protection. Most important of these 
problems are the lack of an extensive data bases and the 
fact that available turbulence expressions, which are empiri- 
cal in nature, have been formulated primarily from overland 
data. 

With the development of systems dependent on laser tech- 
nology, the Department of Defense interest in turbulence in 
this region pertains to assessing and predicting the environ- 
mental influence on optical wave propagation. The refractive 
index problem for the marine boundary layer is more compli- 
cated than the similar problem over land because of the dif- 
ficulties in making point measurements of the fluctuating air 
temperature from which refractive index statistics for optical 
waves are often obtained. For optical waves the atmospheric 
index of refraction variations primary result from tempera- 
ture-induced density changes. These variations cause scintil- 


lation and phase fluctuations of laser beams propagated in the 
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atmosphere and the loss of sight of objects vieved through 
high magnification telescopes. 

An interdisciplinary group at the Naval Postgraduate 
School has been involved in overvater experiments to relate 
the optical properties and the optically relevant meteoro- 
logical parameters of the marine boundary layer to observ- 
able bulk atmospheric parameters. 

The primary purpose of this vork vas to evaluate pre- 
viously formulated overland expressions for the overwater 
regime. Expressions presented by Wyngaard et al. (1971) and 


Friehe (1976) were used to relate C to a stability para- 


T 
meter, either the Richardson number or Z/L where L is Monin- 
Obukhov Length. In addition to this, expressions described 
by Busch and Panofsky (1968) were used to relate € to the 
Stability parameters. 

Ten minute time averaged values of wind, temperature, 
and relative humidity measured at four levels (4.2, 6.6, 
11.3, and 17.7 meters) were examined in this study. The 
most recent and the majority of the observations examined 
were made during the months of September and October 1976 off 
the coast of San Diego, California. Data previously evalua- 
ted by Hughes (1976) and Karach (1976) vere incorporated into 
the final evaluations to supplement the data base. The latter 
data were also measured at four levels (4.2, 6.6, 7.6, and 


ie. Je meters) = 
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e NENA BACKOROUND 


E GENERAL 


C and &€ are optically relevant turbulence parameters be- 


1 
cause they define the intensity of fluctuations in the refrac- 
tive index and the size or scale of refractive index inhomo- 
geneities, respectively. The scales of the temperature 
inhomogeneities that distort image resolution and degrade 
laser beam propagation in the atmospheric surface layer lie 
in the inertial subrange where fluctuations are isotropic. 

The latter characteristic allows theories by Tatarski 
(1961) describing electromagnetic propagation to be applied 


in analyses. Tatarski considered the refractive-index-struc- 


mae parameter, 


9 


which is related to the temperature structure function para- 


meter, — as 


2o -6 P .-2 
Cr 一 (79 د‎ = = 


2 


where P is the pressure and T the temperature of the atmos- 


phere. En may be evaluated in a number of methods. 
The dissipation rate € defines the micro-scale, £o, aS 
Y 
E = Re C2) 
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where y is the kinematic molecular viscosity. The micro- 
scale appears in optical parameters for image resolution. 
For example, the exp(-2) folding distance of the mutual 
coherence, p, is defined as 
_ 1 
Be c 6 py ao 


(1.90. Le K) 


where K is the average wave number. 
2 


T 


The temperature structure-function parameter may be 


RE PERATURE STRUCTURE-FUNCTION PARAMETER, C 


evaluated or estimated by several different methods depend- 
ing on the measurement procedures. At separations r of the 
order of inertial subrange scales, the temperature structure- 


un On in a locally isotropic field has the form 


Pa E T(x+r)] > GoD 
ÓN m 


OM can also be related to the one-dimensional tempera- 
ture spectrum (The Fourier transform of the correlation func- 
tion with separations in the streamwise direction) in the 


inertial subrange as 


Sp 0? - .25 C K (6) 


K. is the streamwise component of vavenumber. 
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A third approach is based on the hypothesis of Corrsin 
(1951) that in the inertial subrange the one dimensional tem- 


perature spectrum can be expressed as 


S. (K) = T age 3 KT SIS (7) 
where x is the rate of molecular destruction of A د‎ 
an Ke SI), e is the rate of molecular dissipation of turbu- 


RON and à is hypothesized to be 


lent kinetic energy (in cm* S 
a universal, dimensionless constant; the same in all flows. 
Equations (6) and (7) yield an expression for cê in terms of 


the rates of dissipation of turbulent kinetic energy and 


temperature variance, viz. 


2 = 
Cm = XE 


173 (8) 


where ß is an empirical constant with a value of 3.20. 


MR MICROSCALE X (OR €) 

Both g and x in Equation (8) represent molecular processes 
and forces which act only over very small distances. As pre- 
vlously indicated, the particular scale length of interest 
in propagation analyses is the microscale of turbulence de- 
fined by Equation (6) and is dependent on only one dynamic 
parameter, the dissipation rate of turbulent kinetic energy, 
e, which can be estimated from the one dimensional spectral 


density of velocity as 


et 219 x7 5/3 (9) 


Su? = € 


where E is an empirical constant determined to be 0.51. 
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PANA BOR SS HYPOTHESIS 

Because measurements are normally made at a point a rela- 
tion between temporal (f) and wavenumber or spatial frequency 
MO às required. According to Lin (1953), if one considers 
that eddies of wavenumber K are correlated only over distances 
not large compared to 27/k, then over such distances the non- 
uniformity of convection velocity should be small relative 
to the convection velocities. Further, if the mean velocity 
is a good approximation to the convection velocity, then 
U >> U'where U' is the spatial gradient of U. If these cri- 
teria are met the turbulence field may be considered as being 
frozen and swept at the mean wind speed U past a stationary 
instrument. This is known as Taylor's hypothesis and it 
allows frequency spectra to be interpreted as streamwise 
wavenumber spectra. This enables the evaluation of C : from 


T 


Equation (6) using point temperature fluctuation measurements. 


2 


E. RELATING Cr AND € RORSTABILITY PARAMETERS 
Turbulence similarity theory relates both e and x to 
measurable stability parameters so that gus can be estimated 


using Equation (8). The 1968 Air Force Cambridge Research 
Laboratory (AFCRL) boundary layer study yielded information 
on the behavior of both g and x and the combination of dimen- 
sional and physical arguments that make up surface-layer 
similarity theory. The AFCRL results demonstrated that the 
Statistics of both the mean and turbulent fields, when pro- 
perly non-dimensionalized, are universal functions of the 


stability parameter (Z/L). 
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Non-dimensionalized e and x values can be expressed as 





EE fi CAES (10) 
U a 
2—— = f, (Z/L) (11) 
PAU. 
L 
where the scaling velocity is U, = (-u'w')*, the scaling tem- 





š E : I 
Peeature is T, = 7 s the Monin Obukhov stability length 
is L = - TU, /kg w'T! and fj(Z/L) and f,(Z/L) are empirical- 


ly determined functions. 
An empirical expression relating e to boundary fluxes 
arises from direct substitution of Equations (10) and (11) 


into Equation (8) as 


)12( )42/1 ,5 ل 


where f,(2/L) is a combination of f, (2/L) and f,(Z/L). The 
formulation of f,(2/L) was examined by Wyngaard, et al. (1971) 
and their results are shown in Figure 1. They obtained the 


following forms for f,(Z/L) 


2727/23 


f2(Z/L) = 4.9[1-7(2/L)] 9e vA (13) 


f,(Z/L) l. 3[1l O CAL] mS LE (14) 


where negative and positive Z/L values correspond to unstable 


and stable conditions respectively. 
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An empirical expression relating e to boundary fluxes is 


Equation (10), where MG - r can replace £, (Z/L), — 


EEX = é@y(2/L) - 2/L (15) 


U,. 
where dy(2/L) is the non-dimensional wind shear. Overwater 
results by Garrett (1972) on the non-dimensional dissipation 
rate, KEZU ` versus Z/L appear in Figure 2. 
Both Equations (12) and (15) can be rewritten so that 


gradient terms appear in the non-dimensional C : and € terms 


Ak 
if the folloving expressions for T, and U, are used 


Ta 7 $-(Z7L) 3z D 
_ KZ gu 
ur = $4 C27L) Fz (17) 
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Fig. 2. The dimensionless dissipation rate versus Z/L. 
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siedo Equation (16) into Equation (12) and Equation 
(17) into Equation (15) yield such expressions for the non- 


dimensional Ee and € terms respectively 


Cm KË 
TEE RA ا‎ — 
OZ mës E 
3 m7 Zea z Z. = 3 
Me iE t |! =) 
2 oz) 


Eee RELATING Z/L TO RICHARDSON NUMBER 
The above expressions can be interpreted more readily by 
relating Z/L to the Richardson number, 


90 


q (—) 
PË MN 
nd) 
05 
as 
Bye = L (RI) (20) 


where f (Ri) is an empirical function which has been well 
documented for the constant flux layer overland (Businger 
et al., 1971) and 9, is the vertical potential temperature, 
(1 + bla). 

The following relationships between Z/L and Ri were 


determined by Dyer and Hicks (1970) for the unstable case 


and by Webb (1970) for the stable case: 


E. sta Vase ct 
Ber 
er — 
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Here a is an empirically derived constant equal to 5. These 
relationships have been recently verified for the overwater 
regime by Paulson, et al. (1972). 

Eembining Equation (18) and Equation (20) yields an ex- 


pression for ony in terms of gradient measurements, 
3507 * fg (RI) (22) 


ne function f (Ri) vhich has been determined for the overland 


constant flux layer is shown in Figure 3 (Wyngaard, et al., 


I). 
Monee oo tr ied 
| Ch 
i X 
° ° pi e À 
10 ° 
= 
` 
@ 
= 
3 Brom wrngeerd, inumi, cad Collins (1971) 
x Out 
sq» 
uU 
Qo Se. 1-.-1 ON icm 
“Le =2O Zi a - 4 '? -iO -08 -06 -04 - 0 2 0 02 04 
Ri 
Fig. 3. The dimensionless temperature-structure parameter 


parameter versus the Richardson number. 
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Similarly, combining Equation (21) and Equation (19) 


yields the following gradient relationships for e 











ë Ni 2 TE š z. A/U 。 
rm ` ROIS R IS =RI(I-ILSR1) I NC (23) 
02 
E 2 Ri -2 
>27, 80,3 ° Ri) 
02 
RI Ri - 3 
I) 
K1 D (24) 


The curve defined by these relationships will be presented 


later (Figure 22). 


G. BULK AERODYNAMIC FORMULA 

An alternate approach for indirect estimates of the tur- 
bulent fluxes in overvater observational experiments is the 
use of bulk aerodynamic formulae. A bulk aerodynamic formula 
is one which relates the boundary fluxes to the wind speed at 
a level and the temperature and humidity differences between 
that level and the surface. Such formulae require several 
assumptions regarding the stability conditions of the boundary 


layer and the turbulent processes within it. 
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Bulk aerodynamic formulae for momentum — sensible 
heat (wT!) and humidity (wg!) fluxes have the following 


forms, 


2 


7 025) 
ee n 35 E D 

w'8 = C, U(T,-T.) = C;, DAT (26) 
e en _ 8 

wg’ = c. 0(Q,-0,) = Cp DAQ (27) 


where the symbols have been defined on page 11. 
Friehe (1976) summarized considerable amounts of recent 
overwater data to obtain overwater estimates of Cp» Cu and 


C He observed C, and C. to be dependent on the values of 


E H E 

UAT as vell as height z. Friehe used the bulk aerodynamic 
formulae to formulate an alternate expression for the estima- 
tion of Z/L in the atmospheric boundary layer above the sea. 


The formulation consisted of replacing the transfer terms 


(T,, and U,) in the expression for L by bulk parameters 
2 2 








(Us = CU, and Tą) to obtain 
Zo. 2317290 + 09272 — 2202 WAT > 25(mKs +) 
L — 2 = 2 
U U 
(25) 
Z ZAT ZAQ - — 
T = -.37115 + 0.132 =) UAT » 25(mKs ^) 


Emerson (28) leads to a new stability parameter X, defined 


ES X = £5 (140.212 A) -20 « DAT « 25(mKs |) 
U 
(29) 
اير‎ = 2202 A) UAT 2 25(mKs 1) 
U 
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Note that the signs of X and X' are the same as 
mL. 
Friehe defined X and X' for the 10 meter level. However, 


these can be generalized to any level by the approximation 


L (20) 
C„(2) Ca CZ) CE CZ) N Z 2 
= = = ss (30) 
C_(10) E (ON C (10) 2 
D H E Ly (7) 
O 


This was done in this study. 


Ape Tetela Equations (28) and (29) into Equation (12) 


yields for -20 < DAT « 25(mKs”?) 
ao 7213 Ae No An 62) unstable 
(AT)? 32 2 x 1073 (1-0.635X) stable 
for UAT > 2e D 
e E -3 -2/3 (31) 
— = 8.03 x 107°(1+2.60x') 
(AT) 
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e ERMENI 


A. PLATFORM AND LOCATION 

Observations aboard the R/V Acania during the experiment 
in September and October 1976 were made in the area from 39° 
15'N to 34°13'N and 123°18'W to 120°20'W; Figure 4. This 


area provided a good sample for open ocean condition. 


be INSTRUMENTATION 
Ten minute averages of temperature, wind, and relative 
humidity were taken at four levels, on two separate masts 
Spatially separated on the forward deck of the ship, Figure 5. 
1. Recording 
Data logging during the experiments was accomplished 
using the NPS developed MIDAS (Microprogrammable Integrated 
Data Acquisition System). This micro-processor based data 
acquisition system utilizes an Intel 8008 central processor 
to control the sampling, averaging, and recording of mean 
meteorological data. All software is written in PL/M to 
Facilitate the writing of self-documenting program. 
كن‎ cam iNeasumement Systems 
Mean wind measurements were made with a Thornthwaite 
Associates cup anemometer wind profile register system, model 
number 104, Figure 6. In operation, the shaft of a three cup 
anemometer unit serves as the shutter between a light source 


and a photocell for each revolution. These three plastic 
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Fig. 5. Mounting arrangements aboard the Acania. 
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Fig. 6. Thornthwaite Associates Cup anemometer 
wind profile register system. 


30 











cone cups are reinforced vith aluminum frames and are mounted 
at 120” intervals on stainless steel tubes attached to the 
main shaft. A characteristic of these cup anemometers is low 
starting speeds vith a small amount of internal friction vhich 
aids in checking inertial overshoot. 

Mean temperature measurements were made with tempera- 
ture sensitive quartz crystal probes, (Hewlett Packard model 
HP-2850) at each level. (Figure 7) RF signals from the 
crystal probes and from a reference oscillator were mixed in 
the HP-2801A readout unit to produce a beat frequency whose 
signature can be analyzed to within 0.001 degrees centigrade 
per hertz. Each sensor simultaneously received pre-experiment 
calibration against a platinum resistance wire thermometer 
in a temperature controlled circulating water bath over the 
expected temperature range. The accuracy in achieving a 0.005 
degree centigrade correction factor was a constant for each 
probe. A 3.7 meter flexible coaxial cable is permanently 
attached to the sensor head and the mast mounted probes are 
housed in an aspirated shelter as seen in Figure 8. Tempera- 
ture values were automatically recorded on a printed tape. 

Mean relative humidity was obtained using the Dunmore- 
lithium chloride sensor and Hygrodynamics Digital I hygrometer 
indicator. (Figure 9) The sensor was also placed in the 
aspirated shelter. The basis of operations of this system is 
resistance charge within an electrolytic solution generating 
a reference voltage variance proportional to the relative 
humidity change. Automatic temperature compensation in the 


Instruments meet the following specifications 
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MW Sr j i ve humidity below 903 
+ 4% relative humidity below 90% 
Calibration of sensors was accomplished by a comparative 
method using a saturated saline solution in a closed con- 
tainer. Observation and recording was accomplished in like 
manner as the temperature. 
3. Fluctuation Measurements 

The temperature fluctuations vere measured using a 
bridge developed at GTE Sylvania, the GTE Sylvania Model 140. 
It was necessary to modify the system slightly for use 
aboard the R/V Acania. 

This system is a small temperature difference sensing 
device originally designed to be carried aloft by a balloon 
and equipped to process and transmit, by pulse-rate modula- 
tion on a radio frequency carrier, this information to a 
ground-base receiving station for demodulation and recording. 

The baseband portion of this system is a balanced 
Wheatstone bridge excited by a 3KHZ signal with a synchronous 
detector on the output. Segments of a very small diameter 
platinum wire serve as temperature sensors in opposite arms 
of the bridge. In the single wire mode, one wire is replaced 
by fixed resistance. An output from the bridge which is 
proportional to the temperature difference between the two 
probes or a temperature change of a single wire in the single 
wire mode is produced. 

The sensor wire is 0.5 cm long and 2.5 x 107°m in 


diameter. This extremely small mass allows a response to 
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temperature variations of up to 1 KHZ while electronic ampli- 
fication allows temperature difference as small as 0.004 
degrees C to be observed, 

Wind fluctuations were recorded on a Honeywell model 
5200 14 channel tape recorder, Real time readout on an eight 
channel chart recorder, Brush model 240, was used to check 
the quality of the signals coming from the sensors. 

tw Tluctuation Statistics 

Spectral analyses to obtain values were performed on 
fluctuating velocity data obtained with Hot Film sensors 
oriented in the vertical, normal to the mean wind. These 
data represent a point measurement and yield temporal or 
time descriptions of the fluctuations. Ten minute segments 
of data originally recorded on magnetic tape were recorded 
into an EMR-Schlumber model 1510 digital spectrum analyzer. 
Procedures for converting spectral values, obtained with the 
analog spectral analyzer, to engineering units and for obtain- 
ing turbulence parameters from the spectra are described in 
the following paragraphs: 

a. oGaling Spectral Plots 

A necessary procedure was scaling the spectral 
plots to relate RMS input voltages to power spectral densi- 
ties (PSD); variance per unit frequency. To obtain power 
spectral density levels, corresponding to RMS voltage inputs, 
calibrated scale charts had to be constructed. 
For purposes of the x-y plot format of the 


analyzer output, the RMS voltages were converted to y = 108,0 
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(voltage) units and a graduated scale was constructed so 
that the logarithm of volts RMS could be interpolated from 
spectral plots. The value of the vertical scale (y) was 
adjusted for each spectrum as a function of both input gain 
and spectral gain. These values were then converted to PSD 
levels for use in calculating € values, 


ir j 3 on between Volts RMS to PSD units used 


was 
ee eos = (cal. level VRMS)* i (32) 
units = 4.5 Bandwidth 
where Bandwidth = Freq. Range 


No. Channels 


à 256 Ep o 
ا‎ 10Hz) 


and cal. level VRMS voltage at y - 0 


na 
II 


ISO put setting). 


Amplitude scaling calibrations vere accomplished 
using externally generated "white noise" signals of 1 Volt 
RMS with a frequency range from 0.1 Hz to 1000 Hz (giving a 
PSD of 107 “V“/Hz). Setting a 3.16 V (10dB) input on the EMR 
1510 Digital Spectrum analyzer insures that 1 Volt corre- 
sponds to y = 0. An example of such a calibration plot for 
a frequency range of 1-200 Hz (the frequency range of inter- 
est) is shown in Figure 10. 


b. Computation of Turbulence Parameters 
from Scaled Spectra 


Values of the turbulence parameter £ were obtained 
from the velocity variance spectra on the basis of the formula 
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for the inertial subrange in vave number space S(k), Equa- 
tion (6), which predicts a -5/3 slope for the spectra when 
plotted in a log-log format. Figure 11 is a typical spectra 
considered in the analyses. Velocity variance spectra con- 
sistently exhibited slopes of -5/3. 

Assuming -5/3 slopes for the variance spectra, 
the intercept of the best fit -5/3 slope 55 the 1 Hz fre- 
quency line was the spectral density denoted (PSD) value used 
in computing the parameter of interest, g. The measured PSD 
value was converted to engineering spectral density units by 


E` relation 


C — ESD 
u 


CA i (33) 


SÉ) 


where C, was the calibration factor for the constant tempera- 
ture, anemometer. Determination of the value for C, has 
been described by Lund (1975). 
Since the velocity fluctuations were measured at 
a fixed point in the flow, the resultant spectral values are 
defined at "temporal" frequencies denoted as S CE) in Equa- 
peon (33). 
5. Analysis of Paired Temperature Sensor 

* values vere estimated by using paired platinum 
wire sensors for variance analysis. This analysis was based 
on the expression for the temperature structure function, 
Equation (5). The analysis is one in which temperature differ- 


ences between the paired resistance wires separated by a 
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distance r is subjected to variance analyses as indicated by 
the term (T'(x) - T'!(x+r))° . 

Voltages corresponding to the temperature differences, 
AT, as measured by the paired sensors vere recorded on analog 
magnetic tape. The variance or root mean square of the tem- 


perature differences required to compute C x was obtained by 


ji 
analyzing the analog voltages with a mean square voltmeter 
(Thermo Systems Incorporated, model 1060) using a 30 second 
time constant for the averaging. The output of the mean-square 
voltmeter (as well as a record of the original signal) were 
recorded on standard strip charts. 

Ten minute averaged values of the RMS voltages were 
abstracted from the strip chart and these values were utilized 
in a program for the IBM 360 computer. This program con- 
verted the mean square voltages values to temperature differ- 


é 2 
ence variance values, from which C were computed. 


E 
ER Pp luctuating Sensor Calibration 


Platinum wire resistance is related to temperature by 


O ee ee (34) 


NN 
where R - resistance, T - calibration temperature, To oo 
Ro = resistance at To: and a = .0036/°C (temperature coeffi- 
cient of resistance for platinum). Each wire's resistance 


was measured in the laboratory at a temperature T, to deter- 
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Temperature fluctuations can be related to vire re- 


sistance fluctuations by differentiating Equation (34) 


dR = aR ar 
or 


1 - 1 
Re aR T (35) 


Recorded voltage fluctuations are then related to re- 
sistance changes by the transfer function, H, in the expres- 


Son 


VS = RS (56) 


where H was determined for the bridges by varying the resis- 
tance of a TSI 1056 Variable Decade Probe Resistance and 
recording the corresponding voltage changes. Voltage vs 
resistance is plotted and the slope of this plot determines 
the transfer function, H. For this experiment H = 10 volts/ 
ohms. 

Finally, to relate recorded voltage fluctuations to 
temperature fluctuations, Equation (36) is substituted into 


Baka tronyg(35) giving 


m M (19:72) 
or 


le En CE (38) 


where Cy 1s the temperature calibration factor in “C0 /volt. 
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Velocity sensors calibrations were accomplished in a 
E ape Using the Hot Wire Calibrator TSI 1125. In this 
process, a known wind velocity 1s blown across the wire and 
the corresponding voltage output is recorded. The sensor wind 


speed iS given by 
u (39) 


where V is the hot-film voltage output and U is the mean wind 
speed for any given level. The constants, a and b, are the 
calibration curve slope and intercept respectively. The con- 
stants, a and b, were derived in the lab. 

To convert the voltage PSD levels to velocity units 


requires a calibration factor given by 
UV (40) 


where C is the calibration factor in m/s/volt, V' is the 
voltage fluctuation, and U' the velocity fluctuation. 


Differentiating Equation (39) yields 


NX” 


EN V (41) 


which when substituted into Equation (40) yields 


1 
2 


NER. = (42) 


The slope of the plotted curve yields the value for a in the 


above equation. 
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GAN ANAL YSIS PROCEDURES 

The data were edited for gross errors or inconsistencies 
due to known instrumentation malfunctions or limitations of 
instruments. An example of limitation on instruments is the 
relative wind direction being greater than 30” off the bow of 
the ship. Earlier experiments have shown that when these con- 
ditions exist, the winds are greatly influenced by the ship 
and therefore are of little value to the experiment. 

The criteria for evaluating or discarding data from in- 
dividual levels or for entire interval depended on performance 
check results obtained during measurements, obvious inconsis- 
tencies, and the availability of fluctuation data. Since mean 
wind U, mean temperature T, and mean humidity a, are expected 
to vary logarithmically with height, these parameters were 
plotted for each time period using the IBM 360 and NPS FORTRAN 
program PLOTP. A best fit straight line was drawn to the data 
points. (See Figure 12.) This procedure was very subjective, 
but was intended to further refine the data by detecting more 
obvious discrepancies and determining if the data did indeed 
vary logarithmically with height. Gradient values of the 
parameters used in actual computations were evaluated using a 
least square fit. By means of the subjective procedure, it 
was determined that the humidity sensor at the second level 
was reading too low during most of the experiment. (See Figure 
lla.) This error was corrected for the least square fit by 


assigning the second level the same value as level three. 
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Beggi. Typical profiles for a) specific humidity — 25 Sept 
75 (0105); b) temperature - 25 Sept 75 (0250); and 
e 012015 - 225° Sept 75 (0551). 
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A. Cy” RESULTS 
l. Gradient Description 
A comparison of 368 overwater Ea and profile results 


with the overland empirical expression, Equation (22), appears 
Pi cure 13, The distribution predicted by Equation (22) is 
the solid line. The circles are mean dimensionless tempera- 
ture structure function parameter (DTSFP) results over .50 
Richardson number (Ri) intervals. The Standard deviation 
corresponding to each mean value appears as vertical lines 
on which the number of values defining the mean appears. Sig- 
nificantly, these results represent more than a factor of 3 
(368 versus 119) of those presented previously. (Davidson, 
et al., 1976). However, the mean DTSFP distribution changed 
very little between the two sets. Larger DTSFP values for 
the overwater regime for near neutral and stable conditions 
(Ri > -.5) can be qualitatively related to the wave influence. 
Decreased DTSFP values for unstable conditions (Ri < -.5) are, 
perhaps, associated with aerodynamic smooth properties for 
the sea surface during light wind conditions. 

A curve describing the observed data was obtained 
using a best fit polynomial analysis. The results of this 


analysis were 
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The influence of the waves can be further examined on 
the basis of observed height variations of — Sami larity 
theory predicts that an has a height variation which depends 
on stability which was discussed by Wyngaard et al. in con- 
junction withan examination of Equation (12). Since T, 15 
independent of height in the surface layer, Equation (12) 
ylelds the following general predictions for the variation of 
a with height based on stability. 

av en near neutral conditions, (Z/L = 0), C : de- 


T 
creases with height as aa Equation (12) becomes 


er = e — (46) 
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IS UD EDL. conditions (Z/L « 0, Equation (13)), 


e decreases more rapidly with height than "S 


EEUU IL, ero would decrease with height as 


d QE since Equation (12) becomes 
2723 x 
c^ B mg. (rry M S (47) 
T K 
a 


with T ke (TFT) independent of height. 


Eno Bs cable conditions (4/1 > Do will decrease 


2 


vith height less rapidly than Z since 


CERN T. 本 E] (48) 


and the term in brackets increases vith increas- 
ing height. 
The vave influence can also be seen in height versus 


uf: er * (computed from cue values) results in Figures 14 


and 15 respectively. In Figure 14, the feature which suggests 


-2/3 


the wave influence is the departure from the Z height dis- 


tribution as the surface is approached. In the upper three 


levels ond appears to fit a ELE 


surface is approached the distribution appears to be better 
-4/3 


distribution, but as the 


defined by a Z Meine diSlEIDULION. In Figure 15, a 


similar change is observed. The upper levels appear to fit a 


mn 41 5 دض‎ 011011 for cë but as the surface is approached 


especially at the one meter level, the C 


fit a pe os mente distribution. 


2 
N values appear to 


Further evidence of wave influence on m appears in 


joint probability-conditional mean results of the normalized 
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temperature variance (oj Ta) over vaves in Figure 16. The 
results and the measurements yielding them have been de- 
Esbed by Davadson (1974). The ratio Em is direct- 
ly related to the DTSFP, defined by Wyngaard, et al. (1971). 
In Figure 16, the vertical axis variable (C/U,) defines wave 
ii Tuence on On/T x C is the phase speed of prevailing swell, 
L is the Monin-Obukhov stability length. The results were 
obtained under stable conditions. The diagonal orientation 
of the contours of the number joint occurrences provides 
qualitative evidence that the waves! influence (C/U,) on 
Or/T:; hence the DTSFP, is as significant as the stability 


w uence (Z/L). 





Na a — + O wv O O 


ka 


Fig. 16. Joint probability-conditional mean results of 
Or/T: for C/U; and Z/L dependence. 
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2. Bulk Aerodynamic 

Results obtained from analyses vith bulk aerodynamic 
variables, Equation (31), are compared to the Wyngaard et al. 
prediction curve in Figure 17. The latter is the solid line. 
The circles are mean dimensionless structure function para- 
meter results over 1.00 X intervals. Standard deviations 
corresponding to mean values appear as vertical lines and the 
number above these vertical lines indicates the number of 
values defining each mean. Larger values of DTSFP are again 
attributed to the wave influence as explained earlier, 

The following was evaluated for the bulk aerodynamic 


expressions using a best fit polynomial analysis. 


ano um stable (49) 


TT 7 (AT) 


.082€.323-X) ^? Unstable (50) 


Three points from an analysis by Friehe (19876) appear 


in Figure 17. A significant aspect of Friehe's analysis of 


2 
T 


fluctuation records which contained large numbers of spikes 


available data was the exclusion of C data from temperature 
exhibiting ramplike signatures, as shown in Figure 18. This 
led to a significant reduction in the amount of data used in 
his comparison. Based on observations by Hicks (1972) and 
hypothesis of Schmitt et al. (1976), Friehe attributed the 
occurrence of the spikes to salt nuclei accumulations on the 
wires with subsequential erroneously sensed temperature 
changes by the system as vater vapor is absorbed or is 


evaporated on the salt nuclei. 
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The occurrence and nature of the spikes were also 
examined vith data collected over Lake Michigan by Davidson 
ANS Elle (1972). Figure 19 includes ramp signatures in 
temperature traces observed over Lake Michigan where salt 
accumulation problem would not be a factor. Therefore, the 
ramp signatures is possibly due to factors associated with 
the nature of the shear flow itself. On this basis there 
did not appear to be sufficient evidence to discard data 
with spikes with ramp signatures. lt is not believed that 
this data being included in this study caused the signifi- 
cant difference between the ae Wyngaard curve and the 
observed data curve. 

A contradiction occurs between results in Figures 13 
l7. In Figure 13, given the same Z and = values for 
both regimes, the plotted curves predict a higher value of 


© for the overwater regime for Ri > -.5 and a higher value 


2 
I 


St in the overland regime for Ri > -.5. On the other 


1 
hand in Figure 17, given the same Z and AT values for both 


regimes, the plotted curves predict a higher value of en im 


the overwater regime for X < 1.5 and a higher value of er 
the ovenland regime for Xx > 1.5. Since Ri and X are cor- 
related they should predict basically the same result. How- 
ever, as is seen in Figures 13 and 17, the only ranges where 
they correlate is the near neutral range. 

In order to explain the reason for this contradic- 


tion, values of Ri vere plotted versus values of X, Figure 20. 


The circles show average values of Ri for X intervals of i.- 
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Temperature trace and velocity fluctuation vectors 
and wave heishts associated with microthermals 
observed over Lake Michigan. 
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Fig. 20. X versus Richardson 
number. 
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Since Z/L is correlated with X, a similar relationship 
should be seen between X and Ri as exists between Z/L and Ri. 
The depicted curve is similar to the latter relationship. 
The conclusion is that the expected correlation existed be- 
tween these two stability parameters. 

A second possible reason for the contradiction is the 


sa 4/3 (30/92)? which norma- 


relationship between Z /AT “and I 
lize ma in the DISFP. Mean values of the former versus the 
latter appear in Figure 21. There should be a positive cor- 
relation between the two and this occurs only for large values 


43 (30/92)7. For small values of the latter there is, 


on 1/2 
however, a negative correlation. This result indicates that 
the two parameters are correctly correlated when small tem- 
perature gradients exist, but are not when the gradients be- 
come large. The contradiction is likely to be due to the 


problem of accurately defining the bulk temperature difference, 


AT, 


ER RESULTS 
1. Gradient Description 
Comparison of the observed e results and the curve 

based on Equations (23) and (24) is shown in Figure 22 where 
the predicted curve is the solid line. Averages over Ri in- 
tervals of .25 of observed data appear as small circles. 

The error bars are standard deviations from the mean within 
each interval. The number of observations defining each 
mean value is at the top of the error bars. The significant 


difference between the predicted curve and the observed data 
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Boren. 20 TS expected on the basis of Equation (23). For 
E 20 the equation is no longer valid since terms in the 
expression become infinite. 

An examination of Equations (13), (14), and (15) 
with the knowledge that U and L are independent of height 
in the surface layer yields the following general prediction 
for the variation of e with heleht. 

a, For neutral conditions, (Z/L = 0), e will decrease 


with height as 77% since Equation (15) becomes 


3 
D 


-p 4 (51)‏ = ع 


DAME OP stablemeonditilons ¡MA 0956 wai deesease 


less rapidly with height than — since Equation 
(15) becomes 
u x 
£ = = 2 (o7 7 Ia) (52) 


c. For unstable conditions, (Z/L < 0), e will de- 


crease less rapidly with height than "E since 
Equation (15) becomes 
UMS i sh 
i gii EIS ) -Z/L C53) 


The quantity in brackets increases with increas- 
ing height. This increase is less than for 
stable conditions. So the decrease of e with 
height for unstable conditions is less than that 


sis table conditions. 
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Results in Figures 23 and 24 indicate that under near neutral 


1 


E ndonse tits a Z height distribution quite well. 


Mere appears to be no change in the height distribution as 


2 
Tp ° 


the surface is approached as there vas vith C 
2. Bulk Aerodynamic 

Results obtained from analyses vith bulk aerodynamic 
variables appear in Figure 25. The predicted curve is the 
solid line, while the circles are mean non-dimensional € re- 
sults over X intervals of 1.0. Again, standard deviations 
Corresponding to mean values are vertical lines and the 
number above these lines indicate the number of values defin- 
ing each mean. The correlation between the non-dimensional 


e values versus X, as described for Z/L by Busch and Panofsky 


MES) is quite good. 


Ck ch EDU PST REGIONS ABOVE SURFACE LAYER 


A Cessna was used in this experiment to extend the data 
base beyond the height of the instrumented towers. Height 


delie i bution of C x values as expressed in Equations (46), 


JE 
(47), and (45), have been shown to be valid to height above 
4000 meters (Tsvang, 1969). Figure 26 shows that a aue S 
distribution exists, indicating stable or near neutral condi- 
tions, from lowest level to about 100 meters. Above 100 
meters, there is a rapid increase of C Values. 

Frisch and Ochs (1975) suggested that if the height oi 


the inversion base, Zi, is another relevant scaling parameter, 


then proceeding in a manner similar to that of Wyngaard et al, 
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ER 23. e- Versus ilog 2. 
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Fig. 24. € versus log 2. 
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Be. 26. > height distribution and temperature profile. 
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Equation (47) could be modified with an additional non-dimen- 


sional function, G(ZIZ.). Equation (47) then becomes 


2 4 = 2/3.2rar\4/3,-4/3 Z 4 
Ca 7 [y1 (1/8) GAN TZ Gi) tou) 
3K al 
where G(ZIZ.) = 1 for ZIZ. << 1 in order to have Equation 


(54) reduce to Equation (47) in the surface layer. 

Frisch and Ochs found a significant departure from a DE RO 
law for values of 2/2; sall. This indicates that the at mos- 
phere is no longer in the "free convection" regime, but is 
affected by the lid on the convective layer at Zs. This 
effect would restrict the integral scale of the turbulent 
temperature fluctuations with height compared to the free con- 
vection value. In addition, because of the layer scale 
circulation and mixing, the temperature variance should be 
greater than the free convection variance. These effects 

vill increase the value of w. relative to the free convec- 

2 


tion prediction. Figure 26 shows the rapid increase of Cr 


values as the inversion is approached. 
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V. CONCLUSIONS 


Application of formulated overland expressions in the 
overwater regime showed a significant difference between ob- 
served cha overwater data and overland predictions. This 
difference is attributed to the wave influence for near 
neutral and stable conditions and to the aerodynamic smooth 
properties of the sea surface during light wind periods for 
moderate unstable conditions. The formulation of these over- 
Water curves wlll greatly enhance the Navy's operational 
ability to describe the marine boundary layer and its in- 
fluence on the optical wave propagation. 

In comparing the DTSFP overland and overwater when plotted 
against Ri and then versus X, significant differences were 
noted. Conclusive explanations for this difference cannot 
be given in this study. Therefore, recommend that serious 
consideration be given to an in-depth examination of this 
significant feature. 

Correlation of non-dimensional e values with the stability 
parameter, Ri, shoved agreement between results and predic- 
tions. For Ri > .20, the theory is not valid and therefore 
there was very poor correlation. Any attempt to describe 
the normalized e for Ri > .20 in the future will require an 
expression which will remain valid for all values of Ri. 


E results appear to be described by the same high dis- 


T 
tributions up to levels vell above the surface layer, 
approximately 60$ of the distance to the marine inversion. 
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